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The hamiltonian of an asymmetric diffusion process with injection and ejection of 
particles at the ends of a chain of finite length is known to be relevant to that 
of the spin-i XXZ chain with integrable boundary terms. However, the inclusion 
of boundary sources and sinks of particles breaks the arrow-reversal symmetry 
necessary for solution via the usual Bethe Ansatz approach. Developments in 
solving the model in the absence of arrow-reversal symmetry are discussed. 



1 Introduction 

The spin-i XXZ Heisenberg chain is the canonical example of an integrable 
quantum hamiltonian. The integrability is assured by the triangle or factori- 
sation equations, 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v). (1) 

Among other interpretations, this relation is often described pictorially in 
terms of the scattering of three particles. This notion of integrability has 
been extended to include the presence of a boundaryJroEl with a similar inter- 
pretation in terms of particles reflecting from a wall. The reflection equation 
is 

R 12 (u - v)Kx{u)R 21 {u + v)K 2 (v) = K 2 {v)R X2 (u + »)*ri(u)J2ai(u - v )- ( 2 ) 

Our interest here lies in the fact that the above R and K matrices define 
the Boltzmann weights for exactly solvable two-dimensional lattice models. In 
particular, the reflection equation follows as the condition for two Sklyanin 
double-row transfer matrices (Fig. 1(a)) to commute. The boundary vertex 
weights can be simply written in terms of the X-matrix elements.a 

Consider the concrete example of the six-vertex model and the related 
XXZ spin chain, for which the R and K matrices are given by 



R = 



( a \ 

sin u sin A 
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K _ - fcsin(£ + it) ^sin(2w) . . . 
' ^sin(2w) fcsin(£ — u) ' 
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(a) (b) 

Figure 1: (a) Sklyanin double-row transfer matrices and (b) the open diagonal lattice. 

where a = sin(A — u^-with k, £, ^, v free parameters. This if-matrix was found 
only quite recentlyflEl 

The vertex weights along the right boundary in Fig. 1 (b) follow as 

^> = fcsin(£ + \u), ^> = yusin(u), 

^> = /csin(£ - \u), ^> = vsm(u). 

In general there is a set of parameters fe±, £±, £i±, v± for each boundary. These 
correspond to the two free ends of the related spin chain, with hamiltoniarJj 

N-l 
3=1 

+ (P- a l + P+ n N + C -°\ + C +°N + d - a l + d + n N) > ( 5 ) 

where a = (a x ,<r v ,a z ) are the Pauli matrices, with cr ± = \(o x ±i<r y ). The 
parameters p± (related to £±) control the strength of the z-component of the 
surface fields. The parameters c± and d± (related to v± and /z±) are the terms 
responsible for breaking the familiar "up-down" Z 2 symmetry. 

Despite some intensive effort, the six- vertex model with boundary weights 
defined by the above if -matrix, and equivalently the XXZ chain with the 
above boundary terms, still defy an exact solution, by which I mean that 
both the transfer matrix and the hamiltonian remain to be diagonalised. The 
major obstacle is readily apparent - the boundary terms arising from the non- 
diagonal elements of the if-matrix break the Zi symmetry by the introduction 
of sources and sinks of arrows or particles. The total spin along a row of 
vertical bonds in the vertex model is no longer a conserved quantity. Such 
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a good quantum number is essential input into the Bethe Ansatz method of 
solution. 

2 Special cases 

2.1 Solution for ji = v = 0. 

Consider first the spiii. .chain, which has been solved for the corresponding 
choice of c± — d± = OpU 

JV-l 

H 1 (A,p^,p+) = J2 (vjVj+i + °s°j+i + A +P- a t +P+°N (6) 

The eigenvalues are given by 

n 

E = (N - 1)A + 4^(cosfcj - A) (7) 

3=1 

with the Bethe equations 

A2(L-i) kj = f-kj (P-, A)/_ fc . (p+, A) Sjh, kj) 

where ff.(a,b) = a — b + e lk and the prime denotes I ^ j. Also S(p,q) — 
s(p,q)/s(q,p), where s(p,q) = 1 — 2Ae 19 + e^ p+q \ This hamiltonian has 
/7 g [su(2)]-symmetryD wlien A = —cos A with p_ = — p + = isinA, known 
also as the 'Potts case'.B 

The six- vertex model on the particular open lattice shown in Fig. 1(b) has 
been considered by a number of authors. For this lattice there is a diagonal-to- 
diagonal transfer matrix as indicated. It was solved by means of the co-ordinate 
Bethe Ansatz for boundary weights corresponding to /i = v = OB The Sklyanin 
double-row transfer matrix (Fig n 1(a)) was also diagonalised for fi = v = 
via the algebraic Bethe Ansatz.13 It is worth noting that a similar double- 
row transfer matrix was diagonalised earlier by Baxter^ who used the Bethe 
Ansatz solution to obtain the surface free energy. Later it was shown how to 
pass from the double-row transfer matrix to the diagonal-to-diagonal_transfer 
matrix by means of a special choice of alternating inhomogeneities.txtj The 
more general double-row transfer matrix was also diagonalised. A key point 
is that the diagonal-to-diagonal transfer matrix does not commute - rather 
the underlying integrability lies in the commutation of the Sklyanin transfer 
matrix with alternating inhomogeneities. A number of other models have also 
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been solved on the lattice in Fig. 1(b) with diagonal if-matricesD Solutions 
of the various quantum invariant spin chains are also known. 

I shall not reproduce the solution of the six-vertex model with diagonal 
isf-matrices here. However, it is worth noting that the parametrisation of 
the diagonal boundary weights in (|j), given via the solution of the reflection 
equation, is particularly convenient when substituted into the corresponding 
weights d and e of Owczarek and Baxteru. 

2.2 The case /x = or v = 0. 

Given the impasse on solving the general problem it came as quite a surprise 
when a solution was reported for the hamiltonianlliHl3 

N-x 

H 2 (a, (3, j, 6) = J2 + ^M+i + 44+1 ) + H s (9) 

3=1 

The surface term H s following from the stochastic dynamics of symmetric 
hopping of particles in one dimension, with particle injection and ejection at 
the boundaries, is 

H s = (a + 7 ) (of - 1) + (7 - - 4) + 

— 2a o-j" + 27 (7+ — (a + 7) — (7 - a)o* + 

2/3 a+ + 25al-(P + S)-((3- 5)a z L (10) 

This is clearly_a special case of the general hamiltonian (|5|). Bethe equations 
were obtainedllj similar in form to those given in (||) with A = 1. However, 
checking the operator algebra in RefJlj reveals that the solution is precisely 
that given in (||), with the identification A = 1, p_ = a ±.7, p+ = (3 + 5. 
Moreover, numerical diagonalisation reveals the equivalence^^ 

H 1 (l,a + j,f3 + 5) = H 2 (a,f3,j,5) (11) 

This is a consequence of the eigenvalues of 

tf 1 (A,p_, P+ ) + d_a 1 + + d + a+ (12) 

being independent of the variables d± (and similarly for crj~ and cr^)lll 

Similar considerations apply to the six- vertex model. Certain combinations 
of sources and sinks at the boundary do not change the eigenspectrum. In 
particular, a boundary sink of arbitrary weight may be included on both edges, 
or equivalently a boundary source on both edges. Such behaviour is consistent 
with inversion relational^ in which the off-diagonal terms involve the prefactors 
fX-V- and n+v + . 
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3 Concluding remarks 

Obtaining the solution with the general iVmatrix remains an outstanding 
problem. It may still be that the solution can be obtained via the pair propaga- 
tion through a vertex technique used for the six- vertex model path antiperiodic 
boundary conditions, where the Z2 symmetry is also broken.ta 

The implications of the special solutions with sources and sinks at the 
boundary to the diffusion problems remain to be fully explored. These include 
the isotropic hamiltonian (^) of relevance, to symmetric hopping with particle 
injection and ejection at the boundarytllta (see also the vertex model in RcfJla). 
It is clear that the corresponding anisotropic XXZ chain enjoys a similar 
property, however its precise meaning in terms of the asymmetric hopping of 
particles needs to be clarified. 
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